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ABSTRACT 

The Descartes circle theorem states that if four circles are mutually tangent in the plane, with 
disjoint interiors, then their curvatures (or "bends") h = ^- satisfy the relation (b\ + 62 + 
63 + 64) 2 = 2{b\ + b\ + b\ + 64). We show that similar relations hold involving the centers of 
the four circles in such a configuration, coordinatized as complex numbers, yielding a complex 
Descartes Theorem. These relations have elegant matrix generalizations to the n-dimensional 
case, in each of Euclidean, spherical, and hyperbolic geometries. These include analogues of 
the Descartes circle theorem for spherical and hyperbolic space. 
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1. Introduction 

We call a configuration of four mutually tangent circles, in which no three circles have a common 
tangent, a "Descartes configuration". The possible arrangements [] of such configurations 
appear in Figure [l], where we allow certain "degenerate" arrangements where some of the 
circles are straight lines. Suppose the radii of the circles are ri,r2,r^,r4. The reciprocals of 
these are the curvatures (or "bends") bj = 1/rj. A straight line is assigned infinite radius; 
then the "bend" is zero. 



In 1643 Rene Descartes [12, pp. 45-50], in a letter to Princess Elizabeth of Bohemia, stated 
a relation connecting the four radii. This relation can be written as a quadratic equation 
connecting the four curvatures, namely: 

Theorem 1.1 (Descartes Circle Theorem) For a Descartes configuration of four mutu- 
ally tangent circles, their curvatures satisfy 

E^(X>) 2 - (L1) 



1 An arrangement of concentric circles having a common tangent is not a Descartes configuration. 



(a) (b) 




(c) (d) 
Figure 1: Descartes configurations 

Descartes considered only the configuration (a) in Figure [|. He did not state the result 



in this form, but gave a more complicated relation algebraically equivalent to (1.1), and his 
sketched proof is incomplete. In 1826 Jakob Steiner (2^, pp. 61-63] independently found the 
result and gave a complete proof. Another independent rediscovery with a complete proof was 
given in 1842 by H. Beecroft and is described in Coxeter ||. Many other proofs have been 
discovered (and rediscovered), some of which appear in Pedoe |23|j . 

The Descartes circle theorem applies to all Descartes configurations of types (a)- (d), pro- 
vided we define the curvatures to have appropriate signs, as follows. An oriented circle is a 
circle together with an assigned direction of unit normal vector, which can point inward or 
outward. If it has radius r then its oriented radius is r for an inward pointing normal and — r 
for an outward pointing normal. Its oriented curvature, (or "signed curvature") is £ for an 
inward pointing normal and — - for an outward pointing normal. By convention, the interior 
of an oriented circle is its interior for an inward pointing normal and its exterior for an outward 
pointing normal. We define an oriented Descartes configuration to be a Descartes configuration 
with the circles having orientations which are compatible in the following sense: either (i) the 
interiors of all four oriented circles are disjoint, or (ii) the interiors are disjoint when all orien- 
tations are reversed. Each Descartes configuration has exactly two compatible orientations in 
this sense, one obtained from the other by reversing all orientations With these definitions, 
the Descartes Circle Theorem remains valid for all oriented Descartes configurations, using 
oriented curvatures. 

In 1936 Sir Frederick Soddy (who earned a 1921 Nobel prize for discovering isotopes) 
published in Nature [27] a poem entitled "The Kiss Precise" in which he reported the result 



above and a generalization to three dimensions. The following year Thorold Gossett [13] 
contributed another stanza giving the general n-dimensional result. We extend the definition 
of a Descartes configuration to consist of n + 2 mutually tangent (n — l)-spheres in M. n in which 

2 The inward pointing orientation of a Descartes configuration is the one in which at least two oriented 
curvatures are strictly positive; the outward pointing orientation is one in which at least two curvatures are 
strictly negative. 
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all pairs of tangent (re — l)-spheres have distinct points of tangency, and orientation is done as 
above. 



Theorem 1.2 (Soddy-Gossett Theorem) Given an oriented Descartes configuration in 
W 1 , if we let bj = 1/rj be the oriented curvatures of the n + 2 mutually tangent spheres, then 



n+2 n+2 



3=1 



The case n = 3 of this result already appears in an 1886 paper of Lachlan [21, p. 498] and his 
proof is given in the 1916 book of Coolidge J?], p. 258]. Thus in calling this result the "Soddy- 
Gossett theorem" we are continuing the tradition that theorems are often not named for their 
first discoverers, cf. Stigler Proofs of the n-dimensional theorem appear in Pedoe [23] and 
Coxeter |J . Pedoe observes that this result is actually a theorem f^] of real algebraic geometry, 
rather than of complex algebraic geometry, in dimensions 3 and above. 

In this paper we present some very simple and elegant extensions of these results, which 
involve the centers of the circles. We show that there are relations, similar to ( |1.2|) , involving 
the centers, together with the curvatures, in the combination curvature x center. Furthermore, 
all these relations generalize to arrangements of re + 2 mutually tangent (n — l)-spheres in re- 
dimensional Euclidean, spherical and hyperbolic spaces, and have a matrix formulation. In the 
process we recover spherical and hyperbolic analogues of the Soddy-Gossett Theorem; these 



were first obtained by Mauldon [22| in 1962. There is a vast literature on this subject, spanning 
two centuries, but (so far) we have not found our matrix formulations in the literature. In 
spirit the ideas trace back at least to Wilker p. 390], see the remark at the end of §4. 



2. The Complex Descartes Theorem 

Given any three mutually tangent circles, with curvatures b%, 62, 63, there are exactly two other 
circles that are tangent to each of these; each gives a four-circle Descartes configuration. See 
Figure ^ for the possible arrangements of the resulting five circles; the three initial circles are 
given by dotted lines. 

The curvatures of these two new circles are the roots of the quadratic equation (|Fl]) (treat- 
ing 64 as the variable). Suppose these roots are 64 and b' 4 . Both can be positive, as in Figure 
2(a), or one may be negative as in Figure 2(b). From ([□]) we have 

b 4 + b' 4 = 2(b 1 + b 2 + b 3 ). (2.1) 

Thus, starting from a Descartes configuration, we can select any one of the four circles and 
replace it by the other circle that is tangent to the remaining three; this gives a new Descartes 
configuration. The new curvature can be obtained from the original four by using (|2.1| ). This 
construction can be repeated indefinitely. We arrive at a packing of circles that fills either (i) 
a single circle, as for example in Figure ||, or (ii) a strip between two parallel lines, or (iii) a 

3 The theorem depends on the fact that the number of real circles, simultaneously tangent to each of n + 1 
mutually tangent real circles with distinct tangents, is exactly two. The total number of complex circles with 
this tangency property is two in dimension n = 2 but typically exceeds two in dimensions n > 3. 
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(a) (b) 

Figure 2: Circles Tangent to Three Tangent Circles 



half-plane, or (iv) the whole plane. Such a figure is called an Apollonian packing, in honor of 
Apollonius of Perga, who considered (about 200 BC) the eight circles that are tangent to each 
of three given circles in general position, cf. Kasner and Supnick j20|. An Apollonian packing 
is completely specified by any three mutually tangent circles in it. 

In constructing the Apollonian packing pictured in Figure 3, we started with four circles 
with vector of curvatures (—1, 2, 2, 3). Each circle has been labelled with its curvature; we notice 



that these are all integers. It is clear from (2.1) that once we have a Descartes configuration 



with all curvatures integral, then in this construction all the curvatures in the packing will be 
integers. 

In 1998, one of us, while computing Figure |3|, with the center of the outer circle located at 
the origin, noticed that the centers of all the circles are rational; in fact in this figure, if a circle 
has curvature b and center (x, y) then (it appeared) bx and by are always integers. Following 



this clue, we were led to the following generalization of (Fl), in which the centers are taken to 
be the complex numbers Zj = Xj + iyj . 

Theorem 2.1 (Complex Descartes Theorem) Any Descartes configuration of four mu- 
tually tangent circles, with curvatures bj and centers Zj = Xj + iyj satisfies 

4 



X>^) 2 = ^(E^) 2 - (2-2) 

3=1 3=1 

Notice that the relation ( p. 2D is of the same form as the original Descartes' relation (|1.1|). 
The complex Descartes theorem implies both the Descartes circle theorem (Fl) and a third 
relation 



4 4 



X>(M»o = s(X>)(I>*i)> (2-3) 
j=i j=i j=i 

These results are obtained by replacing Zj by Zj + w in (|2.2| ), where w is an arbitrary complex 
number, and identifying coefficients of powers of w. 
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Figure 3: An Apollonian packing 
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The complex Descartes theorem also implies a relation similar to (2.1) connecting the 
centers of two circles, each of which is tangent to each of three given mutually tangent circles, 
namely: 



M4 + b' 4 z i = 2 (Mi + b 2 z 2 + hz 3 ). 



(2.4) 



Thus in the iterative construction of an Apollonian packing that we described above, both 
the curvatures and the centers of the new circles can be obtained by simple linear operations 
(followed by divisions). This makes it very easy to draw figures such as Figure 3 using the 
computer. 

The relations in the complex Descartes theorem can be expressed in a more elegant form 
using the matrix 



Q 



2 ■- 



h - 2 l4l 4 



-1 
-1 
1 

-1 1 



(2.5) 



in which l n denotes a column of n l's , and Q2 is the coefficient matrix of the Descartes 
quadratic form 

Q2{XI,X 2 ,X3,X4) := x T Q 2 x = (x\ + x\ + x\ + x\) - ^(xi + x 2 + x 3 + x 4 ) 2 . 

The subscript 2 in Q2 refers to the dimensionality of the space we are considering, 

If b = (61, b 2 , 63, 64) 71 denotes the column vector of curvatures, and c = (bizi, b 2 z 2 , b^z^, b^z^) 1 ' , 
then the Descartes theorem asserts that 



b 1 Q 2 b = 0, 

and the complex Descartes theorem asserts that 

c T Q 2 c = 0. 



(2.6) 



(2.7) 



The complex Descartes Theorem does not completely characterize Descartes configurations in 
the Euclidean plane. There is a slightly stronger result which does, namely: 

Theorem 2.2 (Extended Descartes Theorem) Given a configuration of four oriented 
circles with non-zero curvatures (61, b 2 , 63, 64) and centers {(xi,yi) : 1 < i < 4}, let M 
be the 4x3 matrix 



M :-. 



61 61x1 611/1 

b 2 b 2 x 2 b 2 y 2 

63 b^xz b 3 y 3 

64 64X4 64J/4 



(2.1 



Then this configuration is an oriented Descartes configuration if and only if 

"000 
M T Q 2 M =020 
2 



(2.9) 
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If one or two curvatures bi are zero, and the corresponding centers are infinite, then M can be 
defined in such a way that this matrix identity remains true. 

The complex Descartes theorem follows from this result by applying it to the vector c = 
x+iy, where x and y are the second and third columns of M. The extended Descartes Theorem 
gracefully generalizes to n-dimensions, which we turn to next. 



3. Descartes Configurations in n- Dimensional Euclidean Space 

An n-dimensional oriented Descartes configuration consists of n + 2 mutually tangent oriented 
(n— 1)- spheres Si in n-dimensional space M. n , having distinct tangencies, with the orientations 
compatible in the sense that all interiors are disjoint, either with the given orientation or with 
the reversal of all orientation vectors. Here we suppose that n > 2; the one-dimensional case is 
treated in the concluding section. We often regard a hyperplane as a limiting case of a sphere, 
having zero curvature, with orientation given by a unit normal vector. In what follows an 
"oriented sphere" includes the hyperplane case unless otherwise stated. 

The Soddy-Gossett theorem ( |1.2| ) relates the curvatures of such a configuration of mutually 
tangent n-spheres, and can be written 

Q n (b) := b T Q n b = 0, 

where b = (pi, . . . , b n j r 2) T and Q n ( x ) = x T Q re x is the n-dimensional Descartes quadratic form 
whose associated symmetric (n + 2) x (n + 2) matrix Q n is 

Qn := In+2 l?i+2ln+2- (3-1) 

n 

The Soddy-Gossett theorem has a converse. 

Theorem 3.1 (Converse to Soddy-Gosset Theorem) Ifh = (b\, b n+ 2) T is a 
nonzero real column vector that satisfies 

b T Q n b = 0, (3.2) 

then there exists an oriented Descartes configuration whose oriented curvature vector is b. 
Furthermore any two oriented Descartes configurations having the same oriented curvature 
vector are congruent; that is, there is a Euclidean motion taking one to the other. 

A Euclidean motion is one that preserves angles and distances; it includes reflections, which 



reverse orientations. We do not know an easy proof of this result; a proof appears in [16]. 



The geometry of Descartes configurations is encoded in the curvature vector b. If all bi 
are non-zero and Y^j=i^j > 0> then one of the following holds: (i) all of b\, &2i---) &n+2 are 
positive; (ii) n + 1 are positive and one is negative; (iii) n + 1 are positive and one is zero; or (iv) 
n are positive and equal and the other two are zero. These four cases correspond respectively 
to the following configurations of mutually tangent spheres: (i) n + 1 spheres, with another 
in the curvilinear simplex that they enclose; (ii) n + 1 spheres inscribed inside another larger 
sphere; (iii) n + 1 spheres with one hyperplane (the (n + 2)-nd "sphere"), tangent to each of 
them; (iv) n equal spheres with two common parallel tangent planes. 
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Definition 3.1. Given an oriented sphere S in M n , its curvature- center coordinates consist 
of the (n + l)-vector m(S) given by 

m(5) = (b, bxi, bx n ) (3.3) 

in which b is the signed curvature of S (assumed nonzero) and x(5) = x = (xi,X2, ■■■,x n ) is its 
center. For the "degenerate case" of an oriented hyperplane H its curvature- center coordinates 
m(H) are defined to be 

m(S) = (0,h) (3.4) 

where h := (hi, h%, • • • , h n ) is the unit normal vector giving the orientation of the hyperplane. 

Curvature-center coordinates are not quite a global coordinate system, because they do not 
always uniquely specify an oriented sphere. Given m € M n+1 , if its first coordinate a / then 
there exists a unique sphere having m = m(5). But if a = 0, the hyperplane case, there is a 
hyperplane if and only if = 1> anc ^ m that case there is a pencil of hyperplanes having 

the given value m, which differ from each other by a translation. 

Theorem 3.2 (Euclidean Generalized Descartes Theorem) Given a configuration of 
n+2 oriented spheres Si, S2, ■ ■ ■ S n +2 *n (allowing hyperplanes), let M be the (n+2) x (n+1) 
matrix whose j-th row are the curvature- center coordinates m(Sj), of the j-th sphere. If this 
configuration is an oriented Descartes configuration then 



M T Q n M 




2I n 



diag(0, 2, 2,..., 2). (3.5) 



Conversely, any real solution M to this equation is the matrix of a unique oriented Descartes 
configuration. 

The curvature-center coordinate matrix M of an oriented Descartes configuration deter- 
mines it uniquely even if it contains hyperplanes, because the other spheres in the configuration 
give enough information to fix the locations of the hyperplanes. This result contains the Soddy- 
Gossett theorem as its (1,1)- coordinate. We derive the "if" part of this theorem from the next 
result, proved in §5. However the converse part of this theorem seems more difficult, and we 
do not prove it here. A proof appears in [f[6|] . 

We proceed to a further generalization, which extends the (n + 2) x (n + 1) matrix M to 
an (n + 2) x (n + 2) matrix W obtained by adding an additional column. This augmented 
matrix incorporates information about two oriented Descartes configurations, the original one 
and one obtained from it by inversion in the unit sphere, as we now explain. The definition of 
W may seem pulled out of thin air, but in the next two sections we will show that it naturally 
arises from an analogous result in spherical geometry, which is how we discovered it. 

In n-dimensional Euclidean space, the operation of inversion in the unit sphere replaces 
the point x by x/|x| 2 , where |x| 2 = Y^=i x j- Consider a general oriented sphere S with center 
x and oriented radius r. Then inversion in the unit sphere takes S to the sphere S with center 
x = x/(|x| 2 — r 2 ) and oriented radius f = r/(|x| 2 — r 2 ). Note that if |x| 2 > r 2 , S has the same 
orientation as S. In all cases, 

xx , 
- = - 3.6 

r r 



S 



and 



x| 
r 



r. 



(3.7) 



Definition 3.2. Given an oriented sphere S in R n , its augmented curvature- center coordinates 
are the (n + 2)-vector 



w(S) := (b,b,bxi, ...,bx n ) = (6,m), 



(3.8) 



in which 6 = 6(5'), is the curvature of the sphere or hyperplane S obtained by inversion of S 
in the unit sphere, and m are its curvature-center coordinates. For hyperplanes we define 



w(H) := (b, 0,hx, ...,h n ) = (6,m), 



(3.9) 



where b is the oriented curvature of the sphere or hyperplane H obtained by inversion of H in 
the unit sphere. 

Augmented curvature-center coordinates provide a global coordinate system: no two dis- 
tinct oriented spheres have the same coordinates. The only case to resolve is when S is a 
hyperplane, i.e. 6 = 0. The relation (|3.6[ ) shows that (6, bxi, bx n ) are the curvature-center 
coordinates of S, and if 6 7^ 0, this uniquely determines S, and then, by inversion in the unit 
circle, S. In the remaining case 6 = 6 = then S = S is the unique hyperplane passing through 
the origin whose unit normal is given by the remaining coordinates. 

Given a collection (Si, S2, S n+ 2) of n + 2 oriented spheres (possibly hyperplanes) in W 1 , 
the augmented matrix W associated with it is the (n + 2) x (n + 2) matrix whose j-th row is 
given by the augmented curvature-center coordinates w(Sj) of the j-th sphere. 

The action of inversion in the unit sphere has a particularly simple interpretation in aug- 
mented matrix coordinates. If W is the augmented matrix associated to a Descartes con- 
figuration, and if W is the augmented matrix associated to its inversion in the unit sphere, 
then 



W = W 



1 

1 

I n 



(3.10) 



a result which follows from 



Theorem 3.3 (Augmented Euclidean Descartes Theorem) An oriented Descartes con- 
figuration of n + 2 spheres {Sj : 1 < i < n + 2} in W 1 has an augmented matrix W which 
satisfies 



W 7 Q„W 



0-4 
-4 
24 



(3.11) 



Conversely, any real solution W to this matrix equation is the augmented matrix of a unique 
oriented Descartes configuration. 
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The augmented Euclidean Descartes theorem includes as special cases the "if" direction 
of each of the theorems stated so far, and represents our final stage of generalization of the 
Descartes circle theorem in Euclidean space. In particular, the "if" part of the Euclidean gen- 
eralized Descartes theorem is just ( 3.1 1| ) with the first row and column deleted. In the converse 
direction ^| this theorem gives a parametrization of all oriented Descartes configurations, and 
it is a "moduli space" for such configurations given as an affine real-algebraic variety. 

We discovered the augmented Euclidean Descartes theorem in studying analogues of the 
Descartes theorem in non-Euclidean geometries. In the next section we formulate and prove 
such an analogue in spherical geometry; then in §5 we deduce the augmented Euclidean 
Descartes theorem from it. 



4. Spherical Geometry 

The standard model for spherical geometry S n is the unit n-sphere S n embedded in W l+1 as 
the surface 

S n :={y:yl + yl + ... + y 2 n = l} (4.1) 

with the Riemannian metric induced from the Euclidean metric in IR n+1 by restriction. In this 
model, the distance between two points of E> n is simply the angle a between the radii that join 
the origin of M™" 1 " 1 to the representatives of these points on S n . This distance a always satisfies 
< a < 7T. 

A sphere C in this geometry is the locus of points equidistant (at distance a say) from a 
point in S n called its center. The quantity a = a{C) is the spherical radius or angular radius 
of C; it is the angle at the origin of M n+1 between a ray from to the center of C and a 
ray from to any point of C. Note that there are two choices for the center (and the angular 
radius) of a given sphere; these two choices form a pair of antipodal points of S n . The choice 
of a center amounts to orienting the sphere. In this model the interior of a sphere is a spherical 
cap, cut off by the intersection of the sphere S n with a hyperplane in K n+1 , so (by abuse of 
language) we will also call an oriented sphere a spherical cap. 

The two spherical caps determined by a given sphere are called complementary and the sum 
of their angular radii is it. We define the interior of an oriented sphere to contain all points of 
S n on the same side of the hyperplane as the center of the sphere. If we describe a hyperplane 
by a linear form 

n 

F(y) = J2fm-f, (4.2) 

normalized by the requirement 

i=0 

4 In the converse direction the Augmented Euclidean Descartes theorem is not as strong as the converse in 
the Euclidean generalized Descartes theorem, nor does it imply the converse to the Soddy-Gossett theorem; 
these results require separate proofs. 
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this provides an orientation by denning a positive half-space F(y) > 0. The sphere has center 
f := (/o, ft, ■ ■ ■ , f n ) and has positive radius if and only if |/| < 1. The radius a satisfies 
cos a = /, and the interior of the spherical cap it determines is the region where the linear 
form is positive. A spherical cap can be specified either by a pair (f , a) or by the pair (— f , a— ir), 
while (— f, 7r — a) determines the complementary spherical cap. 

A spherical Descartes configuration consists of n + 2 mutually tangent spherical caps on the 
surface of the unit n-sphere, such that either (i) the interiors of all spherical caps are mutually 
disjoint, or (ii) the interiors of all complementary spherical caps are mutually disjoint. 

Theorem 4.1 (Spherical Soddy-Gossett Theorem) Given a spherical Descartes config- 
uration of n + 2 mutually tangent spherical caps Ci on the n- dimensional unit sphere S n 
embedded in M n+1 , with spherical radius ay subtended by the j-th cap, then these spherical 
radii satisfy the relation 

n+2 n+2 

^(cot«i) 2 = -(^cotai) 2 - 2. (4.3) 

i=\ n 1=1 



This theorem was found by Mauldon [22, Theorem 4] In 1962, as part of a more general 
result allowing non-tangent spheres. He also established a converse: to each real solution of 
( |1.3| ) there corresponds some spherical Descartes configuration, and two spherical Descartes 
configuration with the same data in (fO|) are congruent configurations in spherical geometry. 

The spherical Soddy-Gossett theorem is intrinsic, i.e. it depends only on the Riemannian 
metric for spherical geometry, and not on the coordinate system used to describe the manifold. 
However we shall establish it as a special case of a result that does depend on a particular 
choice of coordinate system. If C is a spherical cap with center y = (yo,yi,U2, ■••2/n+l)> an d 
angular radius a, we define its spherical curvature- center coordinates w + (C) to be the row 
vector 

w+(C) := (cot a, J*_, . . . , ^-). (4.4) 

sin a sin a sin a 

No two spherical caps have the same coordinates w+, since a is uniquely determined by the 
first coordinate, and then the yj are uniquely determined using the other coordinates. 

To any configuration of n + 2 caps C\, . . . , C n +2 we associate the (n + 2) x (n + 2) spherical 
curvature- center coordinate matrix W + whose jth row is w+(Cj). 

Theorem 4.2 ( Spherical Generalized Descartes Theorem) Given a configuration of 
n + 2 oriented spherical caps Cj that is a spherical Descartes configuration, then the (n + 2) x 
(n + 2) matrix W + whose j-th row is the spherical curvature-center coordinates of Cj satisfies 



W^Q n W. 



-2 
2 
2I n 



diag(-2,2,2,...,2). (4.5) 



Conversely, any real matrix W + that satisfies this equation is the spherical curvature- center 
coordinate matrix of some spherical Descartes configuration. 
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The (1, l)-entry of the matrix relation Q4.5| ) is the spherical Soddy-Gossett theorem above. 
This theorem has a remarkably simple proof, which is based on two preliminary lemmas. 
Let J n be the (n + 2) x (n + 2) matrix 



-10 
1 
I n 



diag(-l, 1, 1) 



(4.6) 



Lemma 4.3. (i) For any (n + 2)-vector w +; there is a spherical cap C with w + (C) = w + if 
and only if 

w + J n wl = 1. (4.7) 

(ii) The spherical caps C and C are externally tangent if and only if 

w + (C)J n w + (C') T = -1. (4.8) 

Proof. . (i). If w + comes from a spherical cap with center y and angular radius a, then 

-(cos a) 2 + E]= y] l-(cosa) 2 n 
w+J n w + = — J - = — — = 1 



so (|47D holds. 

Conversely, if Q4.7j ) holds, then one recovers a unique a with < a < tt by setting cot a := 
(w + )i, and one then defines a vector y = (jjq, ...,y n+ %) via yj := ^ , noting that sina 7^ 0. 
The equation ( [4.7| ) now implies that |y| 2 = 1, so y lies on the unit sphere, and we have 
determined a spherical cap giving the vector w+. 

(ii). Two spherical caps with centers y,y' with angular radii a, a' are externally tangent 
if and only if the angle between their centers, viewed from the origin in ]R™ +1 is a + a'. Since 
y and y' are unit vectors, this holds if and only if 

y(y') T = cos(a + a'), 

Now 

w + (C)J n w + (C') T = 7 (-cosacosa / + y(y') T ) 

sin a sin a 

and this gives ( |4.8[ ), using cos (a + a') = cos a cos a' — sin a sin a'. □ 

Lemma 4.4. // A, B are symmetric non-singular n x n matrices and WAW T = B, then 
W T B W = A -1 . 

Proof. . The matrix W is non- singular since B is non-singular. Now invert both sides, then 
multiply on the left by W T and on the right by W. □ 
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Proof of the Spherical Generalized Descartes Theorem. If the caps Cj touch exter- 
nally, we have from Lemma <L3 that 



W+J„W^ 



21 



n+2 



ln+2l 



T 

n+2 



2Qn 1 - 



(4.9) 



Then applying Lemma |4.4j (with A = J n and W = W + ) we obtain 

w£Q n W+ = 2J„~ 1 = 2J n . 



(4.10) 



Conversely, ( 4.10) ) implies (fl.S|), by Lemma 4.4. Looking at the diagonal elements of 
W + J„W^, which are all l's, Lemma |4.3| (i) guarantees that the j-th row of W + is a vec- 
tor w_|_(Cj) for some (uniquely determined) spherical cap Cj, and the off-diagonal elements all 
being —1 shows by Lemma |4.3| (ii) that the caps touch externally pairwise, so form a spherical 
Descartes configuration. □ 



Remark. Wilker [33, pp. 388-390] came tantalizingly close to obtaining the spherical gen- 
eralized Descartes theorem. He termed a spherical Descartes configuration a "cluster", and 



introduced spherical curvature-center coordinates. In a remark he noted our Lemma 4.3 and 
stated equation ( |4.9| ). However he did not invert his formula, via Lemma 4.4 and so failed to 
formulate a result in terms of the Descartes quadratic form. 



5. Stereographic Projection and the Augmented Euclidean Descartes 
Theorem 

We derive the augmented Euclidean Descartes theorem from the spherical Generalized Descartes 
theorem, using stereographic projection. The resulting derivation is reversible, so the spher- 
ical Generalized Descartes theorem and the augmented Euclidean Descartes theorem may be 
viewed as equivalent results. 

Consider the unit sphere in given by Y^oVi = Points on this sphere can be 

mapped into the plane yo = by stereographic projection from the "south pole" (—1,0, ...,0), 
see Figure |3]. (The hyperboloid in the figure will be used later.) 

This mapping (y Q , ...,y n ) -» (x 1} ...,x n ) is given by 

• - " J l<j <n. 



3 i + yo' 

The spherical cap C with center (po, ...,p n ) and angular radius a is the intersection of the 
unit sphere with the plane 

n 

3=1 

The sterographic projection of this cap in the hyperplane yo = is the (Euclidean) sphere S 
with center (x±, ■■■,x n ) and radius r, where 

Xj = , 1 < j < n, 

Po + cos a 



13 



4 pole = (-1,0,...,0) 

Figure 4: Stereo-graphic projection-hyperplane, sphere and hyperboloid 
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and 

sin a 
r = . 

Po + cos a 

If the boundary of the cap C contains the south pole, the corresponding sphere S has infinite 
radius, i.e. it is a hyperplane. 



Proof of the Augmented Euclidean Descartes Theorem. The spherical coordinates of 
the spherical cap C are given by the row-vector 



W + (C) = (cot Q, 



Pi 



Pr. 



sin a sin a 



sin a 



-), 



We relate this to the augmented Euclidean coordinate vector w(S') associated with the corre- 
sponding projected sphere S in the plane yo = 0, given by (|Q|). We have xj/r = pj/sina, 
b = 1/r = cot a + , and we find 



cot a 



Po 



sin a 



Thus 



where 



w(5) = (cot a 



Po 
sin a 



■ , cot a 



Po Pi 



Pr, 



sin a sin a 



sin a 



w+(C)G, 



G 



1 1 
-1 1 
J n 



(5.1) 



Suppose we have a configuration of n + 2 spherical caps Ci, ...,C n+ 2 on the unit sphere. 
These stereographically project into a configuration of Euclidean spheres S\ , . . . , S n +2 m the 
equatorial plane yo = 0, and conversely every configuration of Euclidean spheres lifts to a 
configuration of spherical caps. The map sends spherical Descartes configurations to Euclidean 
Descartes configurations. We assemble the corresponding rows w+(C,-) 3 w(Sj) into matrices 
W + and W, respectively. Then 



W = W+G, 



(5.2) 



and, using the Spherical Generalized Descartes Theorem 4.2, we have 



W T Q n W = G T W^Q n W+G = G T diag(-2,2, ...,2)G 







-4 

2I n 



which proves the augmented Euclidean Descartes Theorem 3.3. □ 
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6. Hyperbolic Geometry 



There are many models of hyperbolic space HI™, of which the three most common are the 
(Poincare) unit ball model, the half space model, and the hyperboloid model. (In two 
dimensions we say "unit disk" and "half-plane" for the first two models.) The unit ball and 
half-space models are described in many places, e.g. Beardon § and Berger [||, Chapter 19]. 
The hyperboloid model, which is less well known, but which is in some ways simpler than the 
others, is described in Beardon Section 3.7], Reynolds [24] and Ryan p5|] . The unit ball and 
half-space models are embedded in K n , though with different metrics, while the hyperboloid 
model is embedded in M n+1 , endowed with a Minkowski metric. Here we need only the unit 
ball and hyperboloid models. A sphere in hyperbolic n-space H n is defined as the locus of 
points equidistant (in hyperbolic metric) from some fixed point in HP, the center. 

The unit ball model consists of the points (y%, . . . ,y n ) in W 1 with Yll=i vf < 1> with the 
ideal boundary being Y^j=x = 1- In this model, the hyperbolic metric is 

n 

ds 2 = (dyf + ■ ■ ■ + dy 2 n )/(l - ^yj) 2 - 

3=1 

and the hyperbolic distance between two points y, y' satisfies 

(n n n \ I n n\ 

(i+Es/i)( 1+ E^ 2 )- 4 E^ / (i-E»?)(i-E^ 2 ) -(6-i) 
3=1 3=1 3=1 ) \ 3=1 3=1 ) 

In this model a hyperbolic sphere (of finite radius) is a Euclidean sphere contained strictly 
inside the unit ball; however its hyperbolic center and hyperbolic radius usually differ from the 
Euclidean ones. 

Points in the hyperboloid model are represented in R n+1 as points on the upper sheet H 7 } 
given by uq > of the two sheeted-hyperboloid i?J cut out by the equation 

Up = 1 + u\ H h Un, 

where H± = U i?" with H™ = —H^t. However the Riemannian metric is not that induced 
from the Euclidean metric on R n+1 , but rather is induced from the Minkowski metric 

ds 2 = —duQ + du\ + ... + du n 

on this space, cf. Beardon [||, p. 49]. The formula for the hyperbolic distance d(u, u') in this 
metric is given by 

cosh(d(u, u')) = uqu' — uxu'i — ... — u n u' n , (6-2) 



see Reynolds [24, formula (6.10)]. One can go between the hyperboloid model and the ball 
model by the change of variables 

Uj 

va = — , for 1 < ? < n, 

^ 1 + UQ ' 
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and in the opposite direction by 

2 2 Vj 
n = — - 1 and Uj = 1 < j < n, 

where 

n 

a = i-E4 

From ( |6.2[) we see that in this model a hyperbolic sphere is represented by the intersection 
of with a hyperplane G(u) = 0, where 

n 

G(u) = g u - ^2 9iUi - g, (6.3) 

i=l 

where 5 > 1 and we require G to be normalized by the requirement that 

n 

So = !+!>!> (M 

i=i 

i.e. the point g := (go,gi,- ■■ ,9n) lies on i?™. This intersection is typically a (Euclidean) 
ellipsoid. The center of the sphere is g, and its radius d satisfies coshd = g. As in the spherical 
case, we define the interior of the hyperbolic sphere to be the region on the same side of the 
plane G(u) =0 as the center. 

If we consider a general hyperplane G(u) = which has the normalized form ( |6.4[ ), inter- 
secting the complete hyperboloid H 1 } U iT™, the intersection may be empty, a single point, an 
ellipsoid (or sphere) on either sheet, a paraboloid on either sheet, or a two-sheeted hyperboloid. 
For later use we term these possibilities virtual hyperbolic spheres, except for the empty set or 
a point. We also assign them an orientation given by the sign of the constant term g in the 
associated linear form G(u). (Replacing G(u) by — G(u) reverses the orientation.) Only the 
points in H 1 ^ correspond to real points in W 1 . In the two-dimensional case parabolas on the 
upper sheet correspond to horocycles; in the unit disc model these are (Euclidean) circles that 
are tangent to the bounding circle. They have infinite radius. In the disc model their centers 
are on the bounding circle, and in the hyperboloid model their centers are at infinity. The 
boundary of the disc model corresponds to a circle at infinity in the hyperboloid model. 

An oriented hyperbolic Descartes configuration is any set of n + 2 mutually tangent oriented 
hyperbolic (n— 1) — spheres in HP, having the property that either (i) all interiors of the spheres 
are disjoint, or (ii) the interiors of each pair of spheres intersect in a nonempty open set. In the 
following result we also allow (oriented) Descartes configurations which include those virtual 
hyperbolic spheres of infinite radius which in the ball model correspond to Euclidean spheres 
lying entirely inside the closed ball and tangent to its boundary. The ideal boundary itself 
forms a single limiting (n — l)-sphere in this sense. 

Theorem 6.1 (Hyperbolic Soddy-Gossett Theorem) The oriented hyperbolic radii 
{sj ■ 1 < j < n + 2} of an oriented Descartes configuration of n + 2 spheres in hyperbolic 
space I™ satisfy the relation 

n+2 n+2 

^(cothsj) 2 = -(^cothsj) 2 + 2. (6.5) 
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This result was found by Mauldon ||22|| . The hyperbolic Soddy-Gossett theorem is intrinsic, 
depending only on the hyperbolic metric. We derive it as a special case of a result which does 
depend on a specific coordinate system, namely that given above for the the hyperboloid model. 

If S is a hyperbolic sphere in with center u = (uq, m, U2, Un+l), and hyperbolic 
radius Sj, we define its hyperbolic curvature- center coordinates w_(S) to be the row vector 

u Ui 



w_(S) := ( coth s, 



u,, 



(6.6) 



sinh s sinh s sinh s 

To a configuration of n + 2 hyperbolic spheres Si, ... , S n+ 2 we associate the (n + 2) x (n + 2) 
matrix W_ whose jth row is w_(Sj). 

Theorem 6.2 (Hyperbolic Generalized Descartes Theorem) Given a configuration of 
(n + 2) oriented hyperbolic spheres which is a hyperbolic Descartes configuration, then the 
associated matrix W whose rows are the hyperbolic curvature- center coordinates of the spheres 
satisfies 



W T Q n W_ 




-2 
2I„ 



dtos(2,-2,2,...,2). 



(6.7) 



The converse of Theorem p.2| does not hold, because some matrices W_ satisfying ( p. 7] ) 
do not correspond to hyperbolic Descartes configurations. We can obtain a converse by allow- 
ing "virtual Descartes configurations" that lie on both sheets of the hyperboloid. One simply 
defines a virtual Descartes configuration to be the image on the two-sheeted hyperboloid re- 
sulting from stereographic projection through (— 1, 0, 0, 0) of any Descartes configuration on 
the unit sphere. The resulting hyperbolic coordinate matrix W_ is to be defined by (|6.8j). One 
can define a "virtual (oriented) hyperbolic sphere", and define its oriented radius and center 
using the formulas following ( |6.3D ; the center and oriented radius of some "virtual hyperbolic 
spheres" may then be (non-real) complex numbers, although by definition their hyperbolic 
curvature-center coordinates will be real. 

Theorem |6.2| is readily deducible from the spherical generalized Descartes theorem via 
stereographic projection through the "south pole" (— 1, 0, 0, ...0) in mapping oriented 

Descartes configurations on the upper sheet H + n of the hyperboloid to spherical Descartes 
configurations which lie entirely on the upper hemispheres^™ of the unit sphere. See Figure 
I in §5. 

In applying stereographic projection, the locus of a hyperbolic (n — l)-sphere on the hyper- 
boloid is mapped to the locus of a spherical (re — l)-sphere on the unit re-sphere, and also to 
the locus of a Euclidean {n — l)-sphere in the plane xq = 0. Note however that the hyperbolic 
center, the spherical center of the associated spherical cap and the Euclidean center of the 
Euclidean sphere are typically all distinct in the sense that they usually lie on three different 
lines through the "south pole" (-1,0, ...,0) in E n+1 . 

One can show that, if W_ is defined as above, and if W + are the spherical Descartes 
coordinates associated to the resulting Descartes configuration on the sphere as in §4, then 
they are related by 



W 



W-i 



10 

1 
I n 



(6.8) 
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We omit a proof of this fact, which can be carried out along the lines of §5. Given it, one 
immediately deduces Theorem 6^2 from Theorem 4.2 , plus a converse if "virtual Descartes 
configurations" are included. 



7. Apollonian Packings 



Using stereographic projection we have a recipe to pass between Euclidean, spherical and hy- 
perbolic Descartes configurations. It gives a one-to-one correspondence between configurations 
W, W+, and W_ given by (5.1) and (|6.8D, namely 



W = Wi 



1 1 
-1 1 

I n 



w 



-1 1 
1 1 
I n 



(7.1) 



Here we have extended the definition of Descartes configuration to "virtual Descartes con- 
figurations" in the hyperbolic case to be configurations on both sheets of the two-sheeted 
hyperboloid. This recipe clearly lifts to Apollonian packings. 

Since the spherical and hyperbolic Soddy-Gossett theorems involve quadratic forms, an 
analogue of the relation ( |2.1[ ) holds in spherical and hyperbolic geometry, permitting the easy 
calculation of the "curvatures" cot a (resp. coths) of circles in spherical (resp. hyperbolic) 
packings. There is a notion of "integral Apollonian circle packing" for such "curvatures" which 
makes sense in spherical and hyperbolic geometry. Furthermore, analogues of the relation ( |2.4| ) 
hold in spherical and hyperbolic geometry as well, permitting the easy calculation of the centers 
in spherical (resp. hyperbolic) Apollonian packings. 

Thus the standard Euclidean Apollonian packing pictured in Figure ^, with center at 
the origin, has a corresponding hyperbolic packing obtained by stereographic projection in 
which the bounding outer circle in the packing is the "absolute" in the unit ball model of the 
hyperbolic plane, and in which the coth(r)'s are all integers, but not the same integers as in 
the Euclidean packing, calculated using W_ in Q. See Figure || 

Those circles that are tangent to the bounding circle are known as horocycles, and have infinite 
hyperbolic radius, so the corresponding value of coth(r) is 1. This explains the large number 
of circles assigned the value 1 in Figure ||, namely all those that touch the outer circle. 

Similarly, in the spherical packing associated to the standard Euclidean packing in Figure 
3, the cot a's are all integers, different from both the Euclidean and hyperbolic cases, starting 
from (0, 1, 1, 2). See Figure |. 

One may notice interesting numerical relations among the integers in these three pack- 
ings. Consider a "loxodromic sequence" of spheres as studied in Coxeter M, |il|] , where 
each sphere is produced by reflection in the largest sphere of the preceding Descartes con- 
figuration, For the "curvatures" one obtains for the Euclidean packing the infinite sequence 
E : (—1,2,2,3,15,38,...), for the spherical packing S : (0,1,1,2,8,21,...) and for the hyper- 
bolic packing H : (—1, 1, 1, 1, 7, 17, ...). Note that S + H = E, since this is so for the initial 
values, and each sequence satisfies the same fourth order linear recurrence relation, which is 
x n+1 = 2x n + 2x n _i + 2x n ^ 2 ~ x n -3, by ([□]). 

In the Euclidean case, there are infinitely many different kinds of Apollonian packings 
having integer curvatures for all circles, see [|17j|. The same occurs for both hyperbolic and 
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Figure 6: A spherical Apollonian packing 



21 



spherical Apollonian circle packings. In the hyperbolic case we also include among such integer 
hyperbolic circle packings some packings which are "virtual packings" . Figure || is generated 
from the basic configuration having coth's (—1, 1, 1, 1), and the next simplest hyperbolic case 
is (-2,3,5,6). 

The Apollonian construction works also in higher dimensions, but gives sphere packings 
only in dimensions two and three; in dimensions four and higher we do not get proper pack- 
ings; after several steps the spheres will overlap, see Boyd ||. However "Apollonian sphere 
ensembles" continue to exist in all dimensions as collections of Descartes configurations, see 



m 



There is a considerable amount of mathematics devoted to circle packings; Kenneth Stephen- 



son's [3C] bibliography of circle-packing papers lists over 90 papers since 1990. For further 
relations of Apollonian packings and the relation of integer Apollonian circle packings to the 
integer Lorentz group 0(1, 3, Z), see our series of papers with Ron Graham and Catherine Yan 
ig,||],fl!], and Soderberg f2§. 

8. Conclusion 

We have extended the Descartes circle theorem, well known for n-dimensional Euclidean space, 
to n-dimensional spherical and hyperbolic space. We presented matrix generalizations of the 
Descartes circle theorem which characterize Descartes configurations in all three geometries, 
and which required for their formulation the use of a particular coordinate system in each of 
these geometries. Mauldon [^] generalized the Soddy-Gossett theorem in all three geometries 
to apply to sets of n + 2 equally inclined spheres, as measured by an inclination parameter 
7, with 7 = — 1 for touching spheres; our matrix theorems can be extended to the case of 
arbitrary 7 as well. 

Interestingly, there are one-dimensional analogues of all these theorems. For the Euclidean 
case in one dimension a "circle" consists of two points bounding an interval, and two "circles" 
are tangent if they have one point in common. The one-dimensional Descartes form is 



-1 -1 
Q i: =l3-l 3 lf= -1 -1 
-1 -1 



A one- dimensional Euclidean Descartes configuration consists of two touching intervals, and 
a third "interval" which is the complement of their union, so that the three intervals cover 
the line H.. Call the third "interval" the infinite interval, and its "length" is defined to be the 
negative of the length of its complement, which is the union of the first two intervals. The 
radius is half the "length." The radii r±,r2,rs, of the three intervals then satisfy 

n+r 2 +r 3 = 0, 

which is equivalent to the Descartes relation 

111 2 2 2 

Qi -,-,- = = 0. (8.3) 

r\ r 2 r 3 77 r 2 nr 3 r 2 r 3 

The value of "curvature x center" of the infinite interval is defined as being equal to the 
"curvature x center" of the finite interval obtained by reflection sending x — > -. This describes 
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a positively oriented Descartes configuration; a negatively oriented one is obtained by revers- 
ing all signs. One can now define a 3 x 3 augmented matrix W exactly as in the augmented 
Euclidean Descartes theorem, and one finds that 



W 7 Q,W 











-4 

2 



•4) 



Conversely, every solution W to this equation corresponds to a one-dimensional Descartes 
configuration. The is even a notion of Apollonian packing in dimension n = 1, but it consists 
of a single Descartes configuration! This holds because there is only a single circle tangent 
to a pair of tangent one-dimensional circles. That is, the Descartes equation (|8.3| ) is linear in 
each curvature variable cij = — separately, instead of quadratic, hence the reflection operation 
which generates new circles to add to the Apollonian packing in dimensions n > 2 does not 
exist. Finally, there are one-dimensional spherical and hyperbolic analogues of these results, 
defined via (7.1), taking n = 1. They can be established by stereographic projection. 

The main results in this paper are theorems in inversive geometry, as described in Wilker |33| , 
also Alexander j|] and Schwerdtfeger [26]. Inversive geometry is the geometry that preserves 
spheres and their incidences, which consists of the study of geometric properties preserved by 
the group M6b(n) of conformal transformations of the space 1" = M n U {oo} ~ S n . The set of 
Descartes configurations form a single orbit under the action of the conformal group, and this 
group appears ^ in our results as the (real) automorphism group 



Aut(Q r 



{N: N T Q n N = Q„} 



of the Descartes quadratic form Q n , which is a Lie group isomorphic to 0(n + 1,1), (see 
Wilker (3^, Corollary p. 390 ] for the isomorphism), and the three generalized Descartes 
theorems given here are invariant under the action of Aut{Q n ). Our results prompt the question: 
Is there a "natural" characterization of the global coordinate systems used in these geometries 
which yields the generalized Descartes circle theorems in the matrix form presented here? 



5 The conformal group is isomorphic to a subgroup of index 2 in Aut(Q n ), and we introduced oriented 
Descartes configurations to keep track of the two cosets. 
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